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ABSTRACT 

In the context of the multiple stellar population scenario in globular clusters (GC), helium 
(He) has been proposed as a key element to interpret the observed multiple main sequences (MS), 
subgiant branches (SGB) and red giant branches (RGB), as well as the complex horizontal branch 
(HB) morphology. In particular, second generation stars belonging to the bluer part of the HB, 
are thought to be more He rich (AY=0.03 or more) but also more Na-rich/O-poor than those 
located in the redder part that should have Y equal to the cosmological value. Up to now this 
hypothesis was only partially confirmed in NGC 6752, where stars of the redder zero-age HB 
showed a He content of Y=0.25±0.01, fully compatible with the primordial He content of the 
Universe, and were all Na-poor/O-rich. Here we study hot blue horizontal branch (BHB) stars 
in the GC NGC 6121 (M4) to measure their He plus O/Na content. Our goal is to complete the 
partial results obtained for NGC 6752, focusing our attention on targets located on the bluer 
part of the HB of M4. We observed 6 BHB stars using the UVES@VLT2 spectroscopic facility. 
Spectra of S/N~150 were obtained and the very weak He line at 5875 A measured for all our 
targets. We compared this feature with synthetic spectra to obtain He abundances. In addition 
O, Na, and Fe abundances were estimated. Stars turned out to be all Na-rich and O-poor and to 
have a homogeneous He content with a mean value of Y=0.29±0.01(random)±0.01(systematic), 
which is enhanced by AY~0.04 with respect to the most recent measurements of the primordial 
He content of the Universe (Y~0. 24-=-0. 25). The high He content of blue HB stars in M4 is 
also confirmed by the fact that they are brighter than red HB stars (RHB). Theoretical models 
suggest the BHB stars are He-enhanced by A(Y)=0. 02-^0. 03 with respect to the RHB stars. The 
whole sample of stars has a metallicity of [Fe/H]=-1.06±0.02 (internal error), in agreement with 
other studies available in the literature. This is a rare direct measurement of the (primordial) He 
abundance for stars belonging to the Na-rich/O-poor population of GC stars in a temperature 
regime where the He content is not altered by sedimentation or extreme mixing as suggested 
for the hottest, late helium flash HB stars. Our results support theoretical predictions that the 
Na-rich/O-poor population is also more He-rich than the Na-poor/O-rich generation and that a 
leading contender for the 2 nd parameter is the He abundance. 



Subject headings: globular clusters: individual (NGC 6121) - Stars: abundances 
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1. Introduction 

In the last few years, following the discovery 
of multiple sequences in the color-magnitude dia- 
grams (CMD) of many globular clusters (GCs), 
the debate on their He content has been re- 
newed. In this respect, the most interesting and 
peculiar clusters are u Centauri and NGC 2808, 
where at least 3 main sequences (MS) are present 
(iBedin et al. Il2004lvmanova et al. 1l2007t|Piotto et al. 



2003). The color of these sequences cannot be ex- 



plain ed in terms of a met allicity effect. In uj Cen- 
tauri Piotto et al. showed that the bluest 



MS is more metal rich than the main red popula- 
tion, while in NGC 2808 they all have the same 
iron content (as inferred from abundances in RGB 
stars). The only remaining parameter affecting 
significantly the position of a star in the MS is 
the helium content, and this was the explanation 



properties of the MS stars in both clusters (INorris 


2004; Piotto et al. 20051: iD'Antona et al. 


2005: 


D'Antona & Ventura 


2007: IPiotto et al. 


2007). 



with the bluer MS stars being more He enriched. 
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Fig. 1.— The CMD of M4 with the observed BHB 
stars indicated as open circles. 

On the other hand almost all GCs observed 



1 Based on observations made with ESO telescopes at La 
Silla Paranal Observatory under programme ID 083.B-0083 



to date ( Carretta et al. 20091) show some kind 
of spread in their element content at the level 
of the RGB, the most evident being the spread 
i n Na and O, elemen ts that are anticorrelated 
(jCarrettaet al. Il2010l) . Na and O abundances 
are also (anti)correlated with other light elemen ts, 



such as C,N,Mg, and Al (|Gratton et al. Il2004h . 



The most natural explanation for this phe- 
nomenon is that clusters experienced an extended 
period of star formation, where the younger pop- 
ulations were born from an interstellar medium 
polluted by products of the CNO (and possi- 
bly NeNa, MgAl) cycle co ming from massive 



stars of the former generation (|Caloi fc D'Antona 
201 ll and reference therein). In this picture the 
interstellar medium is affected by an enhancement 
of its He content, together with N, Na, and Al, 
while C and O turn out to be depleted. This hy- 
pothesis can also explain correlations or anticor- 
rclations of light elements present at the level of 
unevolved stars (jGratton et al. I [20011) . This phe- 
nomenon cannot be due to evolutionary effects like 
deep mixing processes that happen on the red gi- 
ant branch (RGB) only after the first dredge-up 
(i.e. in stars brighter than the RGB-bump), since 
it is also present in MS stars. 

Evidence for a direct correlation between the 
He and Na abundances is now accumulating. 



Bragaglia et al. I ()2010al ) found differences in the 
abundances of Na as well as of other elements 
along the differe nt main sequences of N GC 2808 



In another paper iGratton et al. the same 



group found correlations between several ex- 
pected features likely related to altered He abun- 
dances and the Na abundances for stars along 
the RGB for several clusters (the most clear evi- 
dence being again fo r NG C 2808). Very recently , 



Dupree et al~T(|2010h and IPasquini etalH ([20 lib 



observed that the chromospheric He I 10830 A 
line is stronger in Na-rich RGB stars in u> Cen and 
NGC 2808, respectively, than in Na-poor stars. 

Carretta et al. (2007) found that the extension 



of the Na-0 anticorrelation is related to the ex- 
tension of the horizontal branch (HB), suggesting 
that the anticorrelation may be related also to HB 
morphology. 



Long ago lNorris I (|1981l ) , measuring the strength 
of the CN band at 3839 A in M4, was the 
first to speculate on the possible connection be- 
tween light-element and HB morphology in M4. 
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Catelan fc Freitas Pacheco ( 19951) were the first 
to speculate on the connection between super- 
oxygen-poor stars and blue HB morphology in M3 
and M13. 

It has been clear since the 1960s that the HB 
morphology is related not only to the cluster iron 
content (as expected from the models), but also 
to other parameters (the so called second •parame- 
ter problem) which must account for the fact that 
some clusters have an HB extended or extremely 
extended to the blue, while others of the same 
metallicity do not. 

iGratton et al. I (|2010h showed that the 
main candidate to be the second parame- 
ter is age. However they indicate also that 
at least a third parameter is required, and 
that it is most likely He. 

A spread in He abundance can reproduce 
the HB morphology in GCs (as first noticed by 
Roodlll973l) (see al so the extensive analysis by 



Gratton et al. 2010l ). According to this picture, 



both stars with normal and enhanced He content 
end up on the zero-age horizontal branch (ZAHB) 
after the onset of core He burning. However, if 
mass loss in the previous phase is similar, He-rich 
stars should be less massive (and thus warmer) 
than He-normal ones because they burn H faster 
when they are on the main sequence. So He en- 
hancement could provide the mechanism required 
to spread stars from the redder and cooler HB 
(stars with normal composition) to the hotter and 
bl uer part of the HB (He- rich stars) as discussed 



in lD'Antona et al. I (|2002t ). 

However, the most direct evidence, the spectro- 
scopic determination of He abu ndances directly 



for th e HB stars, is still scarce. IVillanova et al. 
(I2009L hereafter Vi09) analyzed a sample of stars 
in NGC 6752 belonging to the HB in the T off 
range between 8000 and 9000 K. As discussed in 
that paper only stars between 8500 and 11500 K 
are good targets for this purpose because they 
are sufficiently hot to show features of He in 
their spectra but cooler than the Grundahl-jump 



dGrundahl et al. I Il999h . the temperature limit 



above which stars show evidence of metal levita- 
tion and He sedimentati on, which alters th e orig- 
inal surface abundances ( Pace at al. 20061) . Vi09 
found for all ZAHB stars of the redder part of the 
blue HB in NGC 6752 a Na-0 content that, when 
compared with the Na-0 anticorrelation found for 
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Fig. 2. — Example of Na (upper panel) and 
O (lower panel) spectral line fitting for the star 
#45025. Values for Na and O abundances of the 
synthetic spectra are indicated. 



the RGB (jCarretta et al. I l2009h . demonstrates 



that they belong to the Na-poor and O-rich (pri- 
mordial) population. The only evolved HB target, 
which likely comes from the bluer HB, belongs to 
the Na-rich, O-poor population. But, most im- 
portant, they show that stars of the redder ZAHB 
are all Helium-normal (Y=0.25±0.01, where Y is 
the fractional mass of He), in agreement with the 
proposed scenario. 



In addition, very recently ([Marino et al. 11201 if ) 
found a direct correlation between Na and O abun- 
dances and colours of the stars along the HB of 
M4. In fact, RHB stars are all Na-poor and O- 
rich, while BHB are all Na-rich and O-poor. Ac- 
cording to this result we expect the BHB to be 
He-rich, and the RHB to be He-normal. 

While all of this data point toward a connec- 
tion between He spread and HB morphology, a 
key piece of information is still missing to finally 
confirm empirically this scenario. We need to ver- 
ify whether ZAHB stars suitable for He measure- 
ment (8500<T off <11500 K) and belonging to the 
bluer HB indeed show a He-enhancement. The 
aim of this paper is to fill this gap by measuring 
the He content in blue HB (hereafter BHB) stars 
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of M4 (N GC 6121). This clust er has been studied 
in detail (jMarino et al. 1 120081 ). It has a bimodal 
HB with the hotter stars at T c g ^9500 K, and 
a bimodal Na-0 anticorrelation, which is possibly 
related to a spread in He abundance. For this rea- 
son it is the ideal target for our purposes. We 
want to verify whether the hot HB stars of M4 are 
He-rich, as well as Na-rich/O-poor. In Section 2 
we describe the observations. In Sec. 3 we discuss 
the determination of the atmospheric parameters 
of our stars and the line list we used. In Sec. 4, 
5, 6, 7, 8, and 9 we present the spectroscopic and 
photometric analysis, compare our findings with 
the literature, and discuss and summarize our re- 
sults. 




T c „ (K) 

Fig. 3. — T e ff as a function of the B-V color (up- 
per panel) and mass as a function of T e ff (lower 
panel) for our sample of stars (open circles) . Filled 
circles are stars of NGC 6752 from Vi09 for com- 
parison. In the upper panel our data are compared 
with Kurucz's synthetic colors for different gravi- 
ties. Continuous lines are synthetic colors for M4 
metallicity, while dashed lines are synthetic col- 
ors for NGC 6752 metallicity. In the lower panel 
our data are compared with the theoretical mass 
(continuous line) for HB stars. The red cross is 
the mean value for our targets. 



2. Observations, data reduction and mem- 
bership 

Observations of stars in the field of view of M4 
were carried out in June-October 2009 in the con- 
text of the ESO Program ID 083.B-0083. In this 
program we observed a sample of 6 BHB stars 
(see Fig[TJ for a total of 10x45 min. exposures, 
selected from B,V photometry obtained by the 
WFI imager at the ES02.2m telescope. 
The selected targets have spectral type AO ( 
(B - V) ~ 0.0, T G ff - 9000 K ). This choice al- 
lows us to have targets showing He features in 
their spectrum, but not affected by levitation 
or sedimentation. Observations were performed 
with the FLAMES-UVES spectrograph at the 
VLT@UT2(Kueyen) telescope. Spectra of the can- 
didate BHB stars were obtained using the 580nm 
setting with 1.0" fibers, and cover the wavelength 
range 4800-6800 A with a mean resolution of 
R=47000. 

Data were reduced using the UVES pipeline 
(Ballcstcr ct al. 2000), where raw data were bias- 
subtracted, flat-field corrected, extracted and 
wavelength calibrated. Echelle orders were flux 
calibrated using the master response curve of the 
instrument. Finally orders were merged to obtain 
a ID spectrum and the spectra of each star sky- 
subtracted and averaged. Each spectrum has a 
mean S/N^150 per resolution element at 5875 A. 
The membership was established by radial veloc- 
ity measurement. We used the fxcor IRAF utility 
to measure radial velocities. This routine cross- 
correlates the observed spectrum with a template 
having known radial velocity. As a template we 
used a synthetic spectrum calculated for a typi- 
cal A0III star (T off =9000, log(g)=3.00, v t =1.00 
km/s, [Fe/H]=-1.5, roughly the same parameters 
as our stars). Spectra were calculated using the 
2.76 version of SPECTRUM, the local thermo- 
dynamical equilibrium spectral synthesis program 
freely distributed by Richard O. GrajH. 
The error in radial velocity - derived from fxcor 
routine - is less than 1 km/s. Finally, for the 
abundance analysis, each spectrum was shifted to 
rest-frame velocity and continuum-normalized. 
Table [1] lists the basic parameters of the se- 
lected stars: the ID, the J2000.0 coordinates 



2 See |http:// www.phys.appstate.edu/spectrum/spectrum.html 
for more details. 
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(RA fe DEC), V mag nitude, B-V and U-V colors 



(jMomanv et al.l 120031 ). heliocentric radial veloc- 



ity (RVh), T ff, log(g), microturbulence velocity 
(v t , for determination of atmospheric parameters 
see Sections 3 and 4). From the measured radial 
velocities we obtained a mean heliocentric radial 
velocity and a dispersion of: 



< RV H >= 70.2 ± 2.6 km/s, er RV = 6.8 ± 1.8 km/s 



litera- 
(2009): 



The mean velocity agrees well with 
ture values ( see e.g. ISommariva et al7 
< RVh >= 70.3 km/s ). The dispersion we found 
is quite high, but also its error is large. It agrees 
within 2cr with the more recent value derived in the 



literature (jPeterson et al. I Il995l 3.5±0.3 km/s). 



Considering the fact that at the position on the 
CMD of the target stars there is no significant 
background contamination, and that their [Fe/H] 
content agrees very well with the mean value for 
the cluster (see Section 4), we conclude that all 
the observed stars are cluster members. 

3. Atmospheric parameters, rotation and 
chemical abundances 

We used the abundances from Fel/II features 
to obtain atmospheric parameters using the equiv- 
alent width (EQW) method. 

None of our stars show evidence for strong rota- 
tion (see Table [TJ). For this reason EQWs are 
obtained from a Gaussian fit to the spectral fea- 
tures. 

We could measure only a small number of Fe lines 
for each spectrum (5 Fel lines and 9 Fell lines) 
due to the limited spectral coverage. However the 
high quality of our spectra (allowing an accurate 
measurement of the EQWs) and the simultaneous 
use of both Fel/II sets of lines allowed us to obtain 
reliable atmospheric parameters, as confirmed by 
the comparison of theoretical and observational 
results (see below). The analysis was performed 



using the 2009 version of MOOG (|Sneden 1119731 ) 
under a LTE approximat ion coupled w ith AT- 
LAS9 atmosphere models ( Kurucz 19921 ). 



T e ff was obtained by eliminating any trend in 
the relation of the abundances obtained from Fc I 
and Fe II lines with respect to the excitation po- 
tential (E.P.), while microturbulence velocity was 



obtained by eliminating any slope of the abun- 
dances obtained from Fel and Fell lines vs. re- 
duced EQWs. log(g) values were estimated from 
the ionization equilibrium of Fel and Fell lines in 
order to have: 

loge(Fel) = loge(Fell) 

where loge(El.)=log(N E i./N H ) + 12. N Bi . and 
Nh are the density of the element and of hydro- 
gen in number of particles per cm 3 . Adopted val- 
ues for the atmospheric parameters are reported 
in Table [TJ 

All the targets, according to their position on the 
CMD, are consistent with being ZAHB objects 
(see Fig.[Tj. 

The typical random error in T e fj and v t can be ob- 
tained by the following procedure. First we calcu- 
lated, for each star, the errors associated with the 
slopes of the best least squares fit in the relations 
between abundance vs. E.P. and reduced EQW. 
The average of the errors corresponds to the typi- 
cal error on the slopes. Then, we selected one star 
representative of the entire sample (#46061). We 
fixed the other parameters and varied first temper- 
ature and then microturbulence until the slope of 
the line that best fits the relation between abun- 
dances and E.P. or reduced EQW becomes equal 
to the respective mean error. The amount of tem- 
perature and microturbulence variation represent 
an estimate of the random errors, that turned out 
to be 50 K and 0.04 km/s respectively. The error 
in gravity was estimated by satisfying the follow- 
ing equation: 

loge(Fel) - Aloge(Fel) = loge(Fell) + Aloge(Fell) 

where Aloge(Fel) and Aloge(Fell) are the errors 
on Fel/II abundance as given by MOOG. In other 
words we took the values loge(Fel) and loge(Fell) 
that satisfied the ionization equilibrium, decreased 
loge(Fel) and increased loge(Fell) by the errors 
given by MOOG, and estimated a new gravity us- 
ing these new Fel/II abundances. The difference 
with the previous value was assumed to be our er- 
ror on gravity, that turned out the be 0.06 dex on 
average. These errors are to be considered as ran- 
dom and internal. Systematic ones were checked 
as described in the next Section. 

For a detailed description of the linelist we used 
for He, O, Na, and Fe and the solar value we 
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Table 1 

Basic parameters of the observed stars. 



ID 


RA(hh : mm : as) 


DEC(° ■.':") 


V(mag) 


B-V(mag) 


U-V(mag) 


RV H (km/s) 


T off( K ) 


log(g)(dcx) 


v t (km/s) 


vsini (km/a) 


45025 


16:23:34.380 


-26:32:36.60 


13.54 


0.35 


0.85 


80.864 


9030 


3.30 


1.57 


9 


46061 


16:23:47.820 


-26:32:06.00 


13.47 


0.32 


0.77 


75.784 


9250 


3.45 


1.42 


3 


47570 


16:23:34.760 


-26:31:24.70 


13.61 


0.32 


0.77 


66.577 


9370 


3.45 


1.02 


7 


49034 


16:23:37.080 


-26:30:44.60 


13.61 


0.31 


0.73 


62.305 


9500 


3.55 


0.90 


7 


49412 


16:23:36.300 


-26:30:34.00 


13.51 


0.34 


0.84 


74.495 


9170 


3.52 


1.50 


5 


50996 


16:23:27.760 


-26:29:49.00 


13.70 


0.33 


0.78 


66.677 


9120 


3.25 


0.95 
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adopted see Vi09. Suffice to say here that 
and Na abundances were obtained by compari- 
son with synthetic spectra from the O triplet at 
6156-6158 A and the Na doublet at 5889-5895 
A (see Fig. [5] for an example). Na and O are 
known to be affected by NL TE. For this re ason, 
we applied the corrections bv lTakedal (|l997t ) and 
Mashon kina et al. interpolated to the at- 

mospheric parameters of our stars. On the other 
hand, as shown by Vi09, He abundances are not 
affected by NLTE, probably because the He line 
at 5875 A, due to the very high E.P., is formed 
in very deep layers of the atmosphere, where de- 
parture from the LTE condition is negligible due 
to the high density of the gas. A further discus- 
sion is required ab out NLTE correc tion for Fe. 
Some authors (i.e. iQiu et al. 1 120011) claim that 



for A type stars like Vega or our targets a NLTE 
correction of +0.3 dex must be applied to Fel 
LTE abundances, while Fell are formed in LTE 
approximation. Slightly smaller non-LTE correc- 
tio ns have been r ecentl y estimated for such stars 
by iMashonkina ( 201ll ). All our tests show in- 
stead that no NLTE correction is required for Fel, 
at least down to log(g)=3.0. First of all, the anal- 
ysis for Vega done in Vi09 assuming LTE gave us 
the same abundance for Fel and Fell within 0.02 
dex, comparable with the r.m.s scatter. This re- 
sult was confirmed by the analysis of the othe r 



targets of Vi09 and by IVillanova et"al~l (|2009bl ) 



In particular gravities of Vi09 were obtained from 
the wings of H B aimer lines, which are formed in 
LTE approximation. In spite of that no appre- 
ciable difference was found in the mean Fel and 
Fell abundances of the stars (A[Fe/H]=0.04±0.05 
dex). 

There is a further effect to take into account. 
Our results indicate that all the stars are He- 
rich (Y^0.29, see next section). In spite of that 



for our analysis we used atmosphere models with 
normal He-content (Y~0.25). The question is 
whether this difference has some impact on the 
final abundances. We answered this by calcu- 
lating with ATLAS9 a new He-enhanced atmo- 
sphere model for the star #46061, considered as 
representative of our sample. Then we recalcu- 
lated the abundances. We found that O, Na, 
and Fe do not change in a significant way (0.01 
dex or less). On the other hand, the He content 
changes by Aloge(He)=+0.03 dex. The reason 
could be that He lines form deep in the atmo- 
sphere where temperature is higher and the UV 
flux, strongly dominated by H and He opacity, is 
greater. This changes the structure of the atmo- 
sphere in the deepest layers and the strength of 
the spectral lines that are formed there. We took 
into account this effect by applying a correction 
Aloge(He)=+0.03 dex (AY=0.015) to the values 
for He obtained assuming He-normal atmosphere 
models. 

Using spectral synthesis we could also 
measure projected rotational velocities for 
each star. For this purpose we assumed 
a combined instrumental+rotational pro- 
file for spectral features. The instrumen- 
tal profile was assumed to be Gaussian with 
FWHM=R/A (where R is the resolution of 
the instrument). Then vsini was varied in 
order to match the observed profile of Na- 
double lines. Results are reported in Ta- 
ble [^_JTlie_t^riical error on vsini is 1-2 km/s 



( Villanova et al. 2009). The results con- 



firm all these stars are relatively slow rota- 
tors 
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Table 2 

Abundances obtained for the target stars. 



El. 


45025 


46061 


47570 


49034 


49412 


50996 


loge(Ho) 


11.00 


10.95 


11.02 


11.07 


10.94 


11.08 


Y 


0.28 


0.26 


0.30 


0.32 


0.26 


0.32 


[O/Fe] 


0.34 


0.36 


0.30 


0.34 


0.27 


0.38 


[0/Fe] NLTE 


0.23 


0.25 


0.18 


0.21 


0.16 


0.26 


[Na/Fc] 


0.82 


0.85 


0.69 


0.66 


0.73 


0.61 


[Na/Fe] NLT E 


0.47 


0.55 


0.39 


0.37 


0.40 


0.28 


[Fc/H] 


-1.07 


-1.06 


-1.04 


-1.04 


-0.99 


-1.13 



' 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 

0.6 - T 



I O 

• -|- 




0.1 0.2 0.3 0.4 0.5 

[O/Fe] 

Fig. 4. — Na-0 abundances found for our HB tar- 
gets (open circles). Filled circles are data for RGB 
stars from Marino et al. (2008). See Sec. 4 for 
more details. 

4. Systematic errors 

An estimation of the systematic errors (or at 
least upper limits) is very important for our pur- 
pose to compare our He abundance with the pri- 
mordial value of the Universe. 

First of all we checked our T e fr scale compar- 
ing the observed colors with synthetic B,V pho- 
tometry, both from KurucJl. For this purpose 
we as sumed an E(B-V)=0.36 value from lHarris 
(1996). for Vi09 this test was not satisfying be- 
cause temperatures from colors were very different 



3 http: / /kurucz. harvard.edu/ 



from the spectroscopic ones. Meanwhile, we have 
further investigated this problem. By compar- 
ing our photometry, based on observations taken 
with the wide-field imager at the ES02.2m tele- 
scope both for NGC 6752 and M4, with Stetson's 
database (jStetson Il2000l) . we found that the blue 
and red parts of our CMDs are affected by photo- 
metric calibration problems, that can reach several 
hundredths of a magnitude. As a consequence we 
corrected our photometry and the result is plot- 
ted in Fig. [3] (upper panel). Empty circles are 
M4 stars, while filled ones are the old NGC 6752 
data corrected for the reddening of the cluster. 
T e ff and dereddened B-V colors were compared 
with Kurucz synthetic photometry for log(g)=3.00 
(blue) and 3.50 (red), which is roughly the grav- 
ity interval our stars cover. Continuous lines have 
[Fe/H]=-1.00, while dashed lines have [Fe/H]=- 
1.50 in order to cover the metal content of the two 
clusters. We note that colors in this temperature 
regime are very dependent on gravity but are al- 
most unaffected by metallicity. The scatter for M4 
stars is a bit higher than that for NGC 6752, due 
to the differential reddening affecting the former 
cluster. 

If we compare T e fj from B-V colors with those 
obtained from spectroscopy for the two clusters, 
we find that for M4 the latter are underestimated 
by about 100 K, while for NGC 6752 they are 
overestimated by about 200 K. We are left with 
a net difference of 100 K. Part or all of this dis- 
crepancy is related to uncertainties in the value 
we adopted for the interstellar reddening, because 
even a small error of 0.01 mag would lead to an 
error of ~ 150 K in the effective temperatures. 
On average, the good agreement between temper- 



4 http: / /www3.cadc-ccda.hia-iha.nrc- 
enrc. gc. ca /community /STETSON /standards / 
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atures from colours and from spectra is then com- 
forting, and we can assume 150 K as un upper 
limit for the systematic error in temperature. 

A better indication concerning the T e g- system- 
atic errors comes from the comparis on between our 
[Fe/HJ (-1.06, see ne xt section) withlMarino et al. 



(2008, -1.07) and iMarino et all (|201lJ, -1.12). 



First of all we note that both datasets were taken 
with the same instrument, so systematic effects on 
abundances due to the spectrograph (as happens 
for example comparing abundances measured with 
UVES and GIRAFFE, ICarretta et al~ll2009l ) can 
be ruled o ut. The difference is only 0.01 dex with 
respect to Marino et al. j2008l) a nd 0.06 dex with 



respect to Marino et al. I ( 201ll) . On the other 



hand, an error as large as 100 K on our spectro- 
scopic T c g would imply an offset of 0.10 dex. We 
can see that i n the worst case comparison with 
Mar ino et al. ] (1201 lh gives 50 K as an estimate 



for our systematic error in temperature. 

Finally M arino et al. I (|201lh checked the reli- 
ability of their atmospheric parameters compar- 
ing their T e g with tho se derived from models by 
iD'Antona et al. 1 1)20021 ) . Their sample of stars dis- 
tributes with a dispersion of ~50 K and has a null 
offset with respect to the line of perfect agreement 
(see their Fig. 2). This implies a negligible sys- 
tematic error on temparature, valid also for our 
data because in the two papers targets are similar 
and the methodology is the same. 

After these tests we assume AT e ff=50 K as the 
systematic error on our T e ff scale but we will con- 
sider also the upper limit AT c ff <150 K in the dis- 
cussion. 

Systematic errors on our gravity scale can arise 
from the fact that we use Fel/II balance to obtain 
log(g). As said in the previous section this can 
introduce a systematic shift if Fel lines are formed 
in NLTE (while Fell lines form in LTE), as sug- 
gested by some authors. Thanks to our results on 
Vega and on other A stars presented in Vi09, we 
think we have shown that Fel lines can be safely 
treated with LTE approximation, and as a conse- 
quence systematic errors on the log(g) scale should 
be negligible. 

Th is statement is further supported bv lYuce et al. 
(|201lh . Here the authors analize a T eff =12045 K, 



log(g)=3.9 dex star using Kurucz's models, as in 
our case. They obtain atmospheric parameters 



from Stromgren photometry. As they say, AT- 
LAS9 model with the parameters T c g= 1204-5 K, 
log(g)—3.9 derived from the Stromgren photome- 
try meets the requirement of same iron abundance 
from all the different kinds of iron lines. In fact, 
they obtain the same abundance for Fel and Fell 
lines. If this is true for a 12000 K star, where 
NLTE effect (if present) should be larger than in 
our colder stars due to the stronger iron overion- 
ization, we can safely assume that LTE works well 
as an approximation for the atmosphere models of 
our targets, and it can be used to obtain unbiased 
gravities. 

However we decided to perform further tests. 
For this purpose we calculated the mass (jj^-) of 
our stars by inverting the canonical equation: 



log(£)=4.]og(%) 



i°s(fe) + J°s(ft) 



where 



Mi, 



-4.74 



2.5 



and 



Mboi = M v + BC = V - (m - M)-\ 



BC 



Bolometric correction (BC) was taken from 
iFlower ( 19961). and dist ance modulus (to — 



M)v from lHarris 



(|l996h . Results and rela- 
tive errors are reported in Fig. [3] (lower panel) 
and compared with the theoretical model from 



Moni Bidin et al. I (|2007l ). The red cross is the 
mean value for our stars. It agrees very well with 
the theoretical value, well within ler. To match 
exactly the theoretical value we should change our 
log(g) scale by 0.03 dex. 

After these tests we conclude that Alog(g)=0.05 
dex is a reasonable value as the systematic error 
for our gravity scale. 

We finally can estimate an upper limit for the 
systematic error on the He abundance. AT c ff=50 
K gives Aloge(He)=0.02 dex, and Alog(g)=0.05 
dex again gives Aloge(He)=0.02. Those two values 
each translate into AY=0.01 and 0.01 respectively. 
The squared sum gives: 

A sys (Y) = 0.01 (the exact value is 0.014) 

which is our systematic error on the He abun- 
dance. If we consider the upper limit AT e g=150 
K instead, we end up with: 



S 



A sys (Y) < 0.03 
Both values will be used in the discussion. 

5. Results of the spectroscopic analysis 

The chemical abundances we obtained are sum- 
marized in Table [21 From our sample the cluster 
turns out to have a Fe content of: 

[Fe/H] = -1.06 ±0.02 

(internal error only) which ag r ees ve ry well with 



(|2008l) ([Fe/H] = 



the results from iMarino et al. 
1.07±0.01). The agreement is slig htly worse but 
within 0.06 dex with respect to Marino et al. I 
(j201l|)([Fe/H]=-1.12). 



In Fig. 2] we compare Na and O abundances 
for our target stars with the Na- anticorrela- 
tion found by IMarino et al. (120081) for a sample 



of M4 RGB stars. IMarino et al. I (|2008l ) identified 



two separated populations in the Na-0 plane, one 
O-rich and Na-poor, the other O-poor and Na- 
rich. We find that all our blue HB stars have a 
Na/O content which is fully compatible with the 
O-poor/Na-rich population. None of our targets 
belongs to the O-rich/Na-poor group. This is an 
indication that the light-element spread (Na and 
O in this case) is a vital clue to the morphology 
of the HB, and suggests that all O-rich/Na-poor 
objects end-up on the redder part of the ZAHB, 
while O-poor/Na-rich stars end-up on the bluer 
part of the ZAHB. We show this statistically by 
calculating the probability to find by chance six 
stars all belonging to the O-poor/Na-rich popu- 
lation, under the hypothesis that the HB mor- 
phology does not depend on the O/Na chemical 



content. Marino et al. 



(|2008l) show that the two 
populations in M4 contain about 50% of the total 
stars each. Under the previous hypothesis, we ex- 
pect to find the blue HB to have an equal mixture 
of the two populations. The probability of finding 
6 stars, all O-poor/Na-rich, as we found for the 
blue HB stars, is less than 2%. Therefore, we can 
conclude, at 98% confidence level, that the HB po- 
sition does depend on the O/Na chemical content, 
or on a related factor, such as the He abundance. 
A further confirmation of this assessment comes 



(|201lh . In that paper we ana- 
lyzed a large sample of stars of the two HBs in M4. 
For a subsample, observed with UVES, we could 



measure both O and Na, while for the remaining 
stars, observed with GIRAFFE and belonging to 
the red HB, only Na. We found that all the blue 
HB objects are O-poor/Na-rich, as in the present 
paper. On the other hand the two red HB stars 
with measured O are O-rich/Na-poor. All the re- 
maining red HB stars for which we collected GI- 
RAFFE spectra have Na that is compatible with 
the O-rich/Na-poor population. Again, the prob- 
ability to find this result by chance is extremely 
low, in this case well below 1%. 

The clear conclusion is that, at least for M4, the 
HB morphology is related to the light element 
content, specifically Na and O. 

But a more important result regards the He 
content. In Fig. [5] we compare the observed spec- 
tra around the He line with 5 synthetic spectra 
with different He content. The stars turn out to 
have a mean He content of loge(He)=11.01±0.02. 
This translates into: 

Y = 0.29 ±0.01 

In Tab. [2] logaritmic He abundances 
loge(He) for single stars were trasformed in 
mass fraction Y value as well. As noted be- 
fore, our best estimation of the systematic error 
on the He measurement is: 

A sys (Y) =0.01 (or 0.014) 

while the upper limit is: 

A sys (Y) < 0.03 

We conclude that the He content of our stars 
is larger than that of the primordial Universe 
(Y^0. 24-^0. 25) with a level of confidence of about 
4<7 if we only consider the internal error. If we 
combine internal and systematic errors using the 
squared sum, the level of confidence is lower but 
still more than 2a. It is less than 2a (between 1.3 
and 1.6cr, depending on if we assume Y=0.24 or 
Y=0.25 for the primordial content) only if we use 
the upper limit for the systematic error. A pos- 
sible point against the sign ificance of our result 
comes from lSweigart I (|1987l ). This paper suggests 
that after the first dredge-up, the Y content of a 
star increases by ^0.015 for Z=0.001 (roughly the 
metallicity of M4). So we should compare our ab- 
solute abundance with Y~0. 24-hO. 25±0.015~0. 26. 



9 



0.94 r 
0.92 



1 .02 
1 

(j 0.98 
C 0.96 
0.94 
0.92 



1 .02 

1 

M 0.98 
' 0.96 
0.9-1 
0.92 



— 



#45025 



A(A) 



#46061 




x(A) 



#47570 




X(i) 



#49034 




X(A) 



1 

B o.ga 

£ 0.96 
0.94 
0.92 




#49412 



MA) 




x(A) 



Fig. 5. — He lines of our stars. We superimposed on each observed spectrum 5 synthetic ones with 
loge(He)=10. 70,10.80,10.90,11.00,11. 10 respectively. 



This would lower our significance to 1.7a (1.0a 
if we use the upper limit for the systematic er- 
ror). However this He enhancement due 
to the first dredge-up is controversial, be- 
cause other mechanisms (e.g. atomic dif- 
fusion, radiative acceleration, and turbu- 
lence) could be at work and playing a role in 
denning the precise He difference between 
MS and HB. 

While not definitively proven, we believe our 
data strongly hints at a He abundance larger than 
the primordial value or than the value expected 
for a He-normal star after the first dredge-up. In 
the following, we will give various arguments that 
definitely support a similar conclusion. 

First, we can compare the present results with 
the He content of HB stars of NGC 6752 found 
by Vi09. In both studies, we used the same 
spectrograph and methodology to estimate He 
abundance. Therefore, in such a comparison any 
systematic error is cancelled out. Vi09 found 
Y=0.25±0.01 for their sample of stars all belong- 
ing to the redder part of the blue HB of that clus- 
ter. Such a value is compatible with the primor- 
dial He abundance of the Universe. The stars of 



the present study are instead located on the bluer- 
part of the HB of M4, and are then expected to 
be He-enhanced. The Y values of the two sam- 
ples differ at the level of 2.8a. We also applied 
a Kolmogorov-Smirnov test to the two datasets. 
This test concludes that the two distributions are 
not compatible, with a confidence lev el of more 
than 90%. We note that in this case the lSweigart 



(|19871 ) result does not affect the comparison be- 
cause BHB targets in both clusters should have 
experienced the same surface He-enhancement af- 
ter the first dredge-up. Therefore we conclude that 
M4 BHB stars are He-enhanced. 

We next compare M4 and NGC 6752 in more 
detail including O and Na. For this purpose we 
plot in Fig. |5] (upper panel) O and Na abun- 
dances from the followin g sources: M 4 RGB 
stars (filled circles) from iMarino et al. 



NGC 6 752 RGB stars (stars) from lCarretta et al .1 
(|2007ri ). M4 HB stars (open circles, this work), 
and NGC 6752 HB s tars (open squares) from 
Villanova et al. I (|2009h . We see that RGB stars 
of the two clusters define a common Na-0 anti- 
correlation fitted by the curve. NGC 6752 stars 
cover a larger range both in O and Na. HB stars 
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follow the same curve, as expected, but they map 
only a limited part of the anticorrelation, due to 
the fact that they are located in a restricted part 
of the HB. According to the results of this pa- 
per, the Na-0 anticorrelation is accompanied by 
a He-0 anticorrelation. This is shown in Fig. [6] 
(lower left panel) . Red crosses are the mean values 
and errorbars for the two groups of stars, while 
the curve is the fit to the Na-0 trend shifted and 
compressed in the y direction in order to match 
the observed points. This curve represents the 
He content that a stars has according to its Na 
abundance. In Fig. |6] (lower right panel) we report 
also the He-Na correlation. Again red crosses are 
the mean values and errorbars for the two groups, 
while the straight line is the fit to the data. This 
fit has the following form: 

loge(He) = +0.14 ± 0.06 • [Na/Fe] + 10.95 ± 0.02 

In order to verify if these stars have also a ho- 
mogeneous He content, wc performed a detailed 
analysis of internal errors for this element. Helium 
was measured by comparing the observed spec- 
trum with 5 synthetic ones, adopting the value 
that minimizes the r.m.s. scatter of the differ- 
ences. The S/N of the real spectrum and the 
strength of the He line introduce an error in the 
final He abundance which can be estimated cal- 
culating the error on the r.m.s. scatter. For our 
data this error corresponds to uncertainties in the 
abundance of 0.06 dex. This value must be added 
to the errors due to the uncertainties on atmo- 
spheric parameters. As discussed before AT e ff =50 
K gives Aloge(He)=0.02 dex, Alog(g)=0.06 gives 
Aloge(He)=0.02 dex, while the error on microtur- 
bulence has no appreciable influence. The final 
uncertainty on the He abundance is given by the 
squared sum of all the individual errors, and the fi- 
nal result is Aloge(He) to t=0.07 dex. If we compare 
this value to the observed dispersion (0.06±0.02 
dex), we can conclude that our HB stars are com- 
patible with having a homogeneous He content. 

As discussed in the introduction, levitation and 
sedimentation are present in HB stars with tem- 
peratures hotter than 11500 K. As our stars are 
cooler, we expect that they are not affected by 
these phenomena. With the purpose of verifying 
this hypothesis we plotted loge(He), [El/Fe] and 
[Fe/H] vs. T c ff in Fig. [7J For each element we 
plotted the best fitting (continuous) and the ±la 



lines (dashed) . All the elements show a trend that 
is flat within the errors. 

Considering also that o ur abundances agree well 
with the RGB values of iMarino et al. I (|2008l ). we 
conclude that none of the elements measured in 
the present paper show evidence of levitation or 
He sedimentation. 

6. An independent photometric He esti- 
mation 



As noted by many authors (e.g. ID'Antona et al 



2002), the brightness of the HB depends on the He 
content. The higher the He content, the brighter 
the HB. As our BHB stars are He enhanced, while 
the reddest HB stars in each cluster are expected 
to be He-normal (Y~0.25, as we found in NGC 
6752), we also expect a difference in V magni- 
tude between the two branches. The difference 
depends on the exact He difference, but for a 
value of the order of AY~0.04 suggested by our 
analysis, it should be of the order of ^0.15 mag 



(jCatelan et al. I l200sih . In order to verify this 



hypothesis we comp ared our photometry wit h zero 



age HB models by ID'Antona et al. I (|2002h . At 



this point the differential reddening affecting the 
cluster is a problem because it blurs out the exact 
HB location. To solve this issue we corrected the 
C MP for the differential reddening as explained 
~ (|2007l ) 



Sarajedini et al. 



Then, in order to lo- 
cate exactly the HB, we built up the Hess diagram 
plotted in Fig. M 

In this figure, the red HB is clearly visible as a 
bump at B-V=0.93, V=13.49, while the blue HB 
is the overdensity covering the range 0.25 <B-V< 
0.45 and 13.4 <V< 13.8. We overplotted on both 
branches the HB lower envelope lines (the contin- 
uous lines) drawn to be 3xAV brighter than the 
faintest HB star, where AV is the typical photo- 
metric error at the level of the HB. Dashed lines in 
the plot are the zero age HB models for the metal- 
licity of the cluster by ID'Antona et al. I (|2002l) 
with Y=0.24 (fainter line) and Y=0.28 (brighter 
line). Models were shifted to the red by E(B- 
V)=0.36, and shifted in V in order to fit the lower 
envelope of the red HB with the model at Y=0.24. 
The model for Y—0.24 does not fit the lower enve- 
lope of the blue HB, which is brighter, as expected 
if it is He-enhanced. Then, we estimated the 
difference in V magnitude between the model at 
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Fig. 6. — Upper panel: Na-0 anticorrelation as defined by RGB and HB stars of M4 and NGC 6752. Lower 
left panel: He-0 anticorrelation as defined by HB stars of M4 and NGC 6752. Lower right panel: He-Na 
correlation as defined by HB stars of M4 and NGC 6752. See text for more details 



Y=0.24 and the blue HB ZAHB. It turned out 
to be ~0.1 mag. This implies a difference in 
He of the two HBs of AY~0.02 accord ing to 
the models of iD'Antona et al. I (l2002h. and 
of A Y~0.03 according to (Catelan et al. I 
2009T ). 

We conclude the photometric test further 
supports our contention that blue HB stars 
in M4 are He-enhanced (Y=0.29) with re- 
spect to the red HB by A(Y)=u.02^0.03. 

7. Comparison with literature 

Although we claim that M4 blue HB stars are 
He-enhanced, several previou s papers foun d no ev- 



idence for this. The first is Behr (2003). Here 
the author obtained chemical abundances (includ- 
ing He) for a sample of blue HB stars around the 



Grundahl-jump for 6 clusters: M13, M15, M3, 
M68, M92, NGC288. In his Fig. 22 the author re- 
ports [He/H] value as a function of log(T e fj), and 
apparently this does not support our result about 
the He-enhancement because all stars cooler than 
the Grundahl-jump are compatible with a normal 
He content ([He/H]~0). However we wish to call 
attention to the following point. For three clusters 
(i.e. M13, M15, M92) the targets cooler than the 
Grundahl-jump belong t o the reddest par t of th e 



HB so, as in the case of IVillanova et al 



f|2009h . 

they are indeed expected to be He-normal. Two of 
the remaining clusters (NGC288 and M68) do not 
have enough points below the Grundahl-jump to 
derive a firm conclusion. We are left with M3. The 
first impression from Fig. 



22 of Behr 



(120031 ) is 

that these stars have a normal He content, but we 
immediately see that the errors are huge (~0.7 dex 
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Fig. 7. — loge(He), [O/Fe], [Na/Fe], and [Fe/H] vs. T e g for our target stars in M4. For each element we 
plotted the best fitting (continuous) and the ±lcr lines (dashed) 



in the single measurement). Secondly, stars are 
located i n the interm ediate region of the HB (see 
Fig. 1 of lBehr 1120031 ). so their He-enhancement, if 
any, is expected to be close t o the p r imord ial value. 
Finally we point out that all lBehr I (|2003) He mea- 
surements appear to be underestimated. For ex- 
ample M3 has a mean He content of [He/H]~-0.5, 
which translates i nto Y~0.10, c learly too low. The 
explication is that iBehr I ([20031 ) was not interested 
in the absolute abundance of He as we are, but 
rather in the trend of He (among the other ele- 
ments) with position along the HB. So he did not 
check possible systematic errors in his methods 
(e.g. the He linelist). For his purposes this was 
not necessary, but it makes a comparison with the 
present paper very problematic. 

f|2009h . 



A second paper is ICatelan et al 



Here 

the authors compare the HB of M3 with theoret- 



ical models having different He content. They as- 
sert and show in their Fig. 2 and 3 that there 
is no evidence of He enhancement because the 
model with Y=0. 25 fits well the HB. However, 
(|2009h says, the fit is good, ex- 
a small deviation of the lower 
in the immediate 



Catelan et al. 



as 

cept perhaps for 
envelope of the blue HB stars 
vicinity of the "knee" from the theoretical (single- 
Y) ZAHB.. The deviation is of -0.05 mag, with 
the stars brighter than the model. This value is 
small, but clearly visible and points toward a He- 
enhancement of blue HB stars in M3 with respect 
to red HB stars. We would have expected an even 
larger deviation for the bluest HB stars of the clus- 
ter which are expected to be even more He-rich, 
but in that region of the H B the models are degen- 
erate. We can assert that ICatelan et al. ] (|2009l V 
instead of contradicting our finding, actually sup- 
ports it. 
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Finally ISalaris et aTl (|2008h fit the HB 
of NGC 1851 with their theoretical models 
in order to obtain some hint about nature 
of the two populat ions of the clu s ters that 
were identified by Mil one et al. ( 2008h as 
a bimodal SGB. They can fit the HB with 
two models. In the first the two popula- 
tions are assumed to have an age differ- 
ence of 1 Gyrs and the same (primordial) 
He, CNO, anf Fe content. In the second 
they are coeval and have different CNO 
(but the same primordial He and the same 
Fe). Apparently there is no room for an 
He-enhancement. Fi rst of all we notice that 
more recent papers ((Villanova et al." 20ld : 



_ip 

Carretta et al.l 2010bh found that stars in 
NGC 1851 have the same CNO content but 
different Fe (^0.06 dex). A difference in Fe 
have an impact on the age difference that 
is ^0.5 Gyrs (assuming the same CNO) . 



So n one of the models by Sala ris et al. 



(|2008h is appropriate to fit the HB. On the 
other hand if we consider in their Fig. 2 
(e.g. lower panel) the line that connect 
the red part with coolest part of the blue 
synthethic HB and project it on the ob- 
served HB (in order to fit the observed red 
HB), we can see that the observed blue 
HB is slightly brighter than this line. Ac- 
cording to any HB model that assume the 



sam e CNO content (including Sala ris et al. 
2008, see their Fig. 1) this is an indication 



that the blue HB is He-enhancend with re- 
spect the red. This fact is confirmed by a 
recent paper ( Gratton et al.~l 120121 ). Here 
the authors obtained a spectroscopic esti- 
mation of Y=+0.29±0.05 for the BHB. The 
error is quite large, but they show photo- 
metrically that the HB can be fitted only 
assuming Y=+0.248 and Y=+0.280 for the 
red and blue HBs respectivelly. 

We conclude this section by noting that 
literature statements against the presence 
of He-rich stars in GCs are not conclusive, 
and that they are refuted by a new inter- 
pretation of the data or by new results. 



8. Discussion 

Correlations and anticorrelations of chemical el- 
ements observed in GCs (i.e. Na vs. O and Mg 
vs. Al) are attributed to contamination by prod- 
ucts of the H-burning process at high tempera - 



ture (jLanger et al. II1993L iPrantzos et al. 1 12007). 



when N is produced at the expense of O and 
C, and proton capture on Ne and Mg produces 
Na and Al (CNO, N eNa, and MgAl cycles). 



Gratton et al. I (|200l[ ) demonstrated that this con- 
tamination is present also at the level of the MS. 
This means that it is not the result of a mixing 
mechanism present when a star leaves the MS, but 
it is rather due to primordial pollution of the in- 
terstellar material from which stars were formed. 
The mixing mechanism along the RGB can have 
an effect, but only as f ar as C and N are concerned 



((Gratton et al. Il2000h . 



Pollution must come from more massive stars. 
GCs must have experienced some c hemical evolu- 
tion a t the beginning of their life (see lGratton et al. 
20041 for extensive references). The main prod- 
uct of H-burning is He and a He enhancement is 
then expected to be present in stars with an en- 
hancement of N, Na, and Al. The main classes 
of candidate polluters a re: fast-rotating mas- 



sive main-sequence stars ( Decres s in et al. 20071) 



intermediate-mass AGB stars (D'Antona et al 



2002), and also massive binary stars rtde Mink et al 
2009). All these channels can potentially pollute 
the existing interstellar material with pr oducts of 



comp lete CNO burning, including He (see lRenzini 
20081 for an extensive review) . 



Recently IVillanova fe Geisler ( 2011 ) showed 
that for M4 the best candidates that can repro- 
duce the abundance pattern observed for RGB 
stars are massive main-sequence stars. 

In the pollution scenario, a hrst generation 
of O-rich and Na-poor stars (relative to a sec- 
ond stellar generation) is formed from primor- 
dial homogeneous material, which must have been 
polluted by previous SN explosions. This gen- 
eration also has a primordial or close to pri- 
mordial He content (Y=0.24-=-0.25). Then the 
most massive stars of this pristine stellar gener- 
ation pollute the interstellar material with prod- 
ucts of the CNO cycle. This material is kept in 
the cluster due to the relatively strong gravita- 
tional held ( D'Ercole et al. 2008) , and it gives 
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rise to a new generation of O-poor and Na-rich 
stars. This population should also have been He- 
enhanced (Y=0.27-0.35, depending on the ma- 
jor po lluter and on the amount of retained ma - 
terial, iD'Antona et al. 1 120021: iBusso et al~l l2007t) . 
Also the abundance of other elements (including 
s-process elements) may differ in stars of the first 
or second generation. 

In the MS phase, He-rich stars evolve more 
rapidly than He-poor ones, so He/Na-rich (and 
O-poor) stars presently at the turn-off or in 
later phases of evolution should be less mas- 
sive than He/Na- poor (and O-ric h) ones . In 
this framework, ID'Antona et al. 1 ( 2002 ) and 



Carretta et al. I (|2007t ) proposed that a spread in 
He, combined with mass loss along the RGB, may 
be the main ingredient to naturally reproduce the 
whole HB morphology in GCs, as discussed in the 
introduction. According to this scenario, primor- 
dial He-0 rich and Na poor stars end up on the 
redder part of the HB, while stars with extreme 
abundance alterations (strong He enhancement, 
O-poor and Na-rich) may end up on the bluer 
part of the HB, if they have experienced enough 
mass loss during the RGB phase. The extensive 
comparison between the distribution of colours 
along the HB a nd th e Na-0 anti-correlation by 
Gratton et aTl (|2010h also supports this view. 
Also age must play a role, because it determines 
the mean mass of a star that reaches the HB 
(jGratton et al. Il2010t iDotter et al. Il2010l) . 

We stress that the HB represents the ideal locus 
to investigate the effects of chemical anomalies in 
GC stars, as the HB is a sort of amplifier of the 
physical conditions that are the consequence of the 
star composition and previous evolution. 

In this scenario, in M4, we expect that the 
progeny of RGB He-normal, O-rich and Na-poor 
stars should reside in the red part of the HB. 
Therefore, it is not surprising t hat th e observed 



red HB stars of iMarino et al. I (|2011f) are all O- 
rich, Na-poor. But they are too cold to have 
their H e content measured directly. M arino et al. 



(|201ll ) show also that blue HB stars are all O- 
poor and Na-rich. Here we reinforce this result 
with higher S/N spectra and more precise mea- 
surements. But we go further, measuring the He 
content of the blue HB. All our stars turned out to 
be He-enhanced (Y=0.29). By comparing pho- 
tometry with HB models, we showed that the 



blue HB is brighter than the red, as ex- 
pected if they have a different He content. 

If we consider both the He and Na/O content 
of our targets, our res ult strongly confirms th e 
hypothesis sugg e sted b y ID'Antona et al. 1 
Carretta et al. I (|2007t ). iGratton et al 



2002), 



2010D . 



and shows that He (coupled with light-element 
spread) is one of the best candidates (together 
with the metal content and age) to explain the HB 
morphology of GCs and thus is a strong candidate 
for the 2 nd paramete r (or for the 3 rd according to 



Gratton et al. Il2010j ). 



9. Summary 

We studied a sample of BHB stars in M4 with 
a temperature in the range 9000-9500 K, with the 
aim of measuring their He and Na/O content. Tar- 
gets were selected in order to be hot enough to 
show the He feature at 5875 A, but cold enough to 
avoid the problem of He sedimentation and metal 
levitation affecting HB stars hotter than 11500 K. 
Thanks to the high resolution and high S/N of our 
spectra, we were able to measure He abundances 
for all our stars. This is only the second, direct 
measurement of He content from high resolution 
spectra of GC stars in this T G ff regime with the 
aim to test the current models of GC formation 
and multiple-populations. 

Our sample of stars turns out to have [Fe/H]=- 
1.06±0.02, value that agrees well with the liter- 
ature, and belong to the O-poor/Na-rich pop- 
ulation of M4 found in the RGB region by 



Marino et al. 



(120081 ). 

Our targets have a mean He content of Y=0.29 
±0.01 (internal error) which is significantly larger 
than the value found for redder BHB stars in 
NGC 6752 (Y=0.25±0.01), using the same obser- 
vational set up and data reduction procedures. 
Our best estimation for the systematic error is 
AY=0.01 (or 0.014). However it does not affect 
the comparison with NGC 6752 or the photomet- 
ric estimation of the red HB He content. We 
compared also the brightness of the red and the 
blue HB of the cluster, finding that the latter is 
~0.1 mag brighter then the former. This result 
is quite strong, and we can estimate an enhance- 
ment in the He content for the stars in the blue 
vs. red HB of AY=0.02^0.03. This is what we 
expect if stellar position on the HB is driven by 
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Fig. 8. — Hess di agram of the HB. Red a nd blue ZAHB (solid curves) were drawn by hand. Dashed are the 
ZAHB models bv lD'Antona et al. I (|2002h with Y=0.24 (the fainter) and Y=0.28 (the brighter). Single stars 
are shown as red crosses. 



its He and metal content. Our combined ev- 
idence strongly suggests that stars within 
the same globular cluster and among differ- 
ent globular clusters have different He con- 
tent, and that the O-poor/Na-rich population is 
also Hc-cnhanccd with respect the O-rich/Na- 
poor population, as suggested by many theoret- 
ical studies. He is thus a leading contender for the 
2 nd parameter. 

Our results are consistent with theoretical stud- 
ies which predict that, for a given metallicity and 
age, the position of a star in the HB is driven by 
its He content, and that the spread of stars along 
the HB must be related to the Na/O spread visible 
at the level of the RGB. 
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